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Abstract: In order to guarantee stability, known results for MPC without additional terminal 
costs or endpoint constraints often require rather large prediction horizons. Still, stable behavior 

of closed loop solutions can often bo observed even for shorter horizons. Here, wc make use 
of the recent observation that stability can be guaranteed for smaller prediction horizons via 
Lyapunov arguments if more than only the first control is implemented. Since such a procedure 
may be harmful in terms of robustness, we derive conditions which allow to increase the rate 
at which state measurements arc used for feedback while maintaining stability and desired 
performance specifications. Our main contribution consists in developing two algorithms based 
on the deduced conditions and a corresponding stability theorem which ensures asymptotic 
stability for the MPC closed loop for significantly shorter prediction horizons. 
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1. INTRODUCTION 

Model predictive control (MPC), sometimes also termed 
receding horizon control (RHC), deals with the problem of 
approximately solving an infinite horizon optimal control 
problem which is, in general, computationally intractable. 
To this end, a solution of the optimal control problem on 
a truncated, and thus finite, horizon is computed. Then, 
the first part of the resulting control is implemented at 
the plant and the finite horizon is shifted forward in time 
which renders this method to be iteratively applicable, see, 
e.g., Magni et al. (2009). 

During the last decade, theory of MPC has grown rather 
mature for both linear and nonlinear systems as shown 
in AUgower and Zheng (2000) and Rawlings and Mayne 
(2009). Additionally, it is used in a variety of industrial 
applications, cf. Badgwell and Qin (2003) due to its capa- 
bility to directly incorporate input and state constraints. 
However, since the original problem is replaced by an 
iteratively solved sequence of control problems - which 
are posed on a finite horizon - stability of the resulting 
closed loop may be lost. In order to ensure stabilty, several 
modifications of the finite horizon control problems such as 
stabilizing terminal constraints or a local Lyapunov func- 
tion as an additional terminal weight have been proposed, 
cf. Keerthi and Gilbert (1988) and Chen and Allgower 
(1999). 

Here, we consider MPC without additional terminal con- 
tstraints or costs, which is especially attractive from a 
computational point of view. Stability for these schemes 
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can be shown via a relaxed Lyapunov inequality, see Griine 
(2009) for a set valued framework and Griine and Pannek 
(2009) for a trajectory based approach. Still, stable behav- 
ior of the MPC controlled closed loop can be observed in 
many examples even if these theoretical results require a 
longer prediction horizon in order to guarantee stability. 
The aim of this paper consists of closing this gap between 
theory and experimental experience. In particular, we de- 
velop an algorithm which enables us to ensure stability for 
a significantly shorter prediction horizon and, as a con- 
sequence, reduces the computational effort in each MPC 
iteration significantly. The key observation in order to ex- 
tend the line of arguments in the stability proofs proposed 
in Griine (2009) is that the relaxed Lyapunov inequality 
holds true for the MPC closed loop for smaller horizons if 
not only the first element but several controls of the open 
loop are implemented, cf. Griine et al. (2010). However, 
this may be harmful in terms of robustness, cf. Magni 
and Scattolini (2007). Utilizing the internal structure of 
consecutive MPC problems in the critical cases along the 
closed loop trajectory, wc present conditions which allow 
us to close the loop more often and maintain stability. 
Indeed, our results are also applicable in order to deduce 
performance estimates of the MPC closed loop. 

The paper is organized as follows: In Section 2, the MPC 
problem under consideration is stated. In the ensuing sec- 
tion, we extend suboptimality results from Griine and Pan- 
nek (2009) and Griine et al. (2010) in order to give a basic 
MPC algorithm which checks these performance bounds. 
In Section 4, we develop the proposed algorithm further 
in order to increase the rate at which state measurements 
are used for feedback and state our main stability theorem. 
In the final Section 5 we draw some conclusions. Instead 



of a separated example, we use a numerical experiment 
throughoiit the paper which shows the differences and 
improvcanents of the presented result in comparison to 
known results. 

2. PROBLEM FORMULATION 

We consider nonlinear discrete time control systems 

x{n+l) = f{x{n),u{n)), x{0) = Xq (1) 

with x{n) e X c X and u{n) e U C J/ for n e No 
where No denotes the natural numbers including zero. The 
state space X and the control value space U are arbitrary 
metric spaces. As a consequence, the following results are 
applicable to discrete time dynamics induced by a sampled 
- finite or infinite dimensional - system, see, e.g., Altmiiller 
et al. (2010) or Ito and Kunisch (2002). Constraints may be 
incorporated by choosing the sets X and U appropriately. 
Furthermore, we denote the space of control sequences 
u:No-^UbyUNo. 

Our goal consists of finding a static state feedback u = 
n{x) € U which stabilizes a given control system of type 
(1) at its unique equilibrium x* . In order to evaluate the 
quality of the obtained control sequence, we define the 
infinite horizon cost functional 

oo 

Joo{xo,u) = '^i{x{n),u{n)) (2) 

n=0 

with continuous stage cost ^ : X x U — >■ with £{x* ,0) — 
and i{x,u) > for x ^ x* . Here, denotes the 
nonnegative real numbers. The optimal value function is 
given by Voo{xo) = inf^gnj^o Jca{x(),u) and, for the sake of 
convenience, it is assumed that the infimum with respect 
to u G U^" is attained. Based on the optimal value function 
Voo{-), an optimal feedback law on the infinite horizon can 
be defined as 

Hoo{x{n)) := argmin{V'oo(a;(n + 1)) + £{x{n),u)} (3) 

using Bellman's principle of optimality. However, since the 
computation of the desired control law (3) requires, in gen- 
eral, the solution of a Hamilton-Jacobi-Bellman equation, 
which is very hard to solve, we use a model predictive 
control (MFC) approach instead. The fundamental idea of 
MFC consists of three steps which are repeated at every 
discrete time instant: 

• Firstly, an optimal control {uq, u^, . . . , u*j^_ _i) G 
for the problem on a finite horizon [0, N), i.e., 

N-l 

Jn{xo, u) = ^ i{xu{k; xo), u{k)) (4) 
fc=o 

is computed given the most recent known state xq 
of the system (1). Here, Xu{-;xo) corresponds to the 
open loop trajectory of the prediction with control u 
and initial state xq- 

• Secondly, the first element /XAr(a;o) = Uq of the 
obtained sequence of open loop control values is 
implemented at the plant. 

• Thirdly, the entire optimal control problem consid- 
ered in the first step is shifted forward in time by 
one discrete time instant which allows for an iterative 
application of this procedure. 



The corresponding closed loop costs are given by 

oo 

Joo{xo,fJ'N) = ^i{x^j,{n),iJ,N{x^^{n))) 

n=0 

where Xfij^{-) denotes the closed loop solution. We use 

Un{-;xo) = argmin JAr(xo, w) (5) 

to abbreviate the minimizing open loop control sequence 
and V/v(.xo) = min^gjjA' Jjv(.to, ?i) for the corresponding 
optimal value function. Furthermore, we say that a con- 
tinuous function p : ]R>o — )• ]R>o is of class /Coo if it 
satisfies p(0) = 0, is strictly increasing and unbounded, 
and a continuous function l3 : M>o x M>o — > M>o is of class 
ICC if it is strictly decreasing in its second argument with 
limt^oo /3(r, t) = for each r > and satisfies /3(-, t) £ /Coo 
for each t>0. 

Note that, different to the infinite horizon problem (2), 
the first step of the MFC problem consists of minimizing 
the truncated cost functional (4) over a finite horizon. 
While commonly endpoint constraints or a Lyapunov 
function type endpoint weight are used to ensure stability 
of the closed loop, see, e.g., Keerthi and Gilbert (1988), 
Chen and AUgowcr (1999), Jadbabaie and Hauser (2005) 
and Graichen and Kugi (2010), we consider the plain 
MFC version without these modifications. According to 
Griine and Fannek (2009), stable behavior of the closed 
loop trajectory can be guaranteed using relaxed dynamic 
programming. 

Proposition 2.1. (i) Consider the feedback law /ijv : X — >■ 
U and the closed loop trajectory a;^„(-) with initial value 
a;(0) = Xo e X. If 

V^r(a;^„(n)) > VN{Xf,^{n + I)) + al{x^,^{n), iiN{x{n))) 

(6) 

holds for some a e (0, 1] and all n e No, then 

olVoo{xq) < aJoo{X0,IJ'N) < Vn{xo) < Voo{xo) (7) 

holds. 

(ii) If, in addition, there exist ai,a2,oi3 e /Coc such that 
Q!i(||a;||) < Vn{x) < a2(||a:|l) and £{x,u) > a:i{\\x\\) hold 
for all X = x^jy(n) G X, n G Nq, then there exists /3 G ICC 
which only depends on ai,Q;2,a3 and a such that the 
inequality ||a;^„(n)|| < /3(||a;o||, n) holds for all n G No, 
i.e., a;^„(-) behaves like a trajectory of an asymptotically 
stable system. 

The key assumption in Froposition 2.1 is the relaxed 
Lyapunov-inequality (6) in which a can be interpreted 
as a lower bound for the rate of convergence. Prom the 
literature, it is well known that this condition is satisfied 
for sufficiently long horizons N, cf. Jadbabaie and Hauser 
(2005), Grimm et al. (2005) or Alamir and Bernard (1995), 
and that a suitable A'^ may be computed via methods 
described in Fannek (2009). To investigate whether this 
condition is restrictive, we consider the following example. 

Example 2.2. Consider the synchronous generator model 

Xl{t) = X2{t) 

X2{t) = -61X3 (t) sin xi(t) - 62X2 (i) + P 

X3{t) = 63 COSXi(t) — &4X3(t) + E + u{t) 

with parameters bi = 34.29, 62 = 0.0, 63 = 0.149, 
64 = 0.3341, P = 28.22, and E = 0.2405 firom Galaz et al. 
(2003) and the initial value xo = (1.02,0.1,1.014). 



For Example 2.2 our aim is to steer the system to its 
equilibrium x* « (1.12,0.0,0.914). To this end, we use 
MPC to generate results for sampling time T = 0.1 and 
running cost 

e{x,u)= [ \\ip{t;x,u) - x*f + X\\u{t)\\'^ dt 
Jo 

with picccwise constant control u{t) = u for t € [0, T) and 
A = 10~^ where i^{-;x,u) denotes the solution operator of 
the differential equation with initial value x and control u. 
For this setting, we obtain stability of the closed loop via 
inequality (6) for N > 30, cf. Figure 1. However, Figure 
1 shows stable behavior of the closed loop for the smaller 
prediction horizon iV = 19 without significant deviations 
- despite the fact that the relaxed Lyapunov-inequality 
from Proposition 2.1 is violated at three instances. Our 




Fig. 1. Comparison of closed loop solutions for iV = 30 
(solid) and A'' = 19 (dashed) where (□) marks a 
violation of the suboptimality bound a = 0.1 

numerical experiment shows that, even if stability and the 

desired performance estimate (7) cannot be guaranteed 
via Proposition 2.1, stable and satisfactory behavior of the 
closed loop may still be observed. Hence, since the compu- 
tational effort grows rapidly with respect to the prediction 
horizon, our goal consists of developing an algorithm which 
allows for significantly reducing the prediction horizon. To 
this end, we adapt the thereotical condition given in (6). 
Moreover, the proposed algorithms check - at runtime - 
whether these conditions are satisfied or not. 

Note that, since we do not impose additional stabilizing 
terminal constraints, feasibility of the NMPC scheme is an 

issue that cannot be neglected. In particular, without these 
stabilizing constraints the closed loop trajectory might run 
into a dead end. To exclude such a scenario, we assume the 
following viability condition to hold. 

Assumption 2.3. For each a; € X there exists a control 
u e U such that f{x, u) e X holds. 

3. THEORETICAL BACKGROUND 

In Section 2 we observed a mismatch between the relaxed 
Lyapunov inequality (6) and the results from Example 
2.2 which indicates that the closed loop exhibits a stable 
and - in terms of performance - satisfactory behavior 
for significantly smaller prediction horizons N . In Griine 



et al. (2010) it has been shown that the suboptimality 
estimate a is increasing (up to a certain point) if more 
than one element of the computed sequence of control 
values is applied. Hence, we vary the number of open 
loop control values to be implemented during runtime at 
each MPC iteration, i.e., the system may stay in open 
loop for more than one sampling period, in order to deal 
with this issue. Doing so, however, may lead to severe 
problems, e.g., in terms of robustness, cf. Magni and 
Scattolini (2007). Hence, the central aim of this work is the 
development of an algorithm which checks - at runtime - 
whether it is necessary to remain in open loop in order to 
guarantee the desired stability behavior of the closed loop 
by suboptimality arguments, or if the loop can be closed 
without loosing these characteristics. 

To this end, we introduce the list S = (s(0), s(l), . . .) C 
No, which we assume to be in ascending order, in order 
to indicate time instances at which the control sequence 
is updated. Moreover, we denote the closed loop solution 
at time instant s(n) by a;„ = Xf^f^{s{n)) and define m„ := 
s{n + 1) — s(n), i.e., the time between two MPC updates. 
Hence, 

^f^N (sin) + rrin) = a;^„ (s(n + 1)) = x„+i 
holds. This enables us - in view of Bellman's principle of 
optimality - to define the closed loop control 

l4f{-,Xn) ■■= a-rgmm\VN-mn{xu{m„;Xn)) (8) 

m„ — 1 

+ ^ i{xuik;xn),u{k))y 

In Griinc ct al. (2010), a suboptimality degree ajv^m^ 
relative to the horizon length N and the number of controls 
to be implemented rn„ has been introduced in order to 
measure the tradeoff between the infinite horizon cost 
induced by the MPC feedback law IJ.%{-; •), i.e. 

OC' m,,, — 1 

V^^ixo) := 51 E ^ (x^sjk;xn),l^%{k;xn)) , (9) 

n=0 k=0 

and the infinite horizon optimal value function V„c{')- In 

particular, it has been shown that given a controllability 
condition, i.e., for each xq there exists a control u^g such 
that 

^{xu„„{n;xo),Ua;o{n;xQ)) < Ca" mm£{xo,u) (10) 

holds with C > 1, cr <E (0, 1), then there exists a prediction 
horizon N G N>2 and to S {1, . . . , [iV/2j} such that 
0(N,m > S for an arbitrarily specified a € (0, 1). 
Example 3.1. Suppose (10) to hold with C = 4 and a = 
0.6 and fix a = 0.275. Then the pair (TV, to) = (15,6) 
imphes aN,fn ~ 0.294 > a, cf. Figure 2, i.e., the closed- 
loop satisfies the performance bound a. However, N = 25 
is the smallest prediction horizon which yields ajv.m > S 
for m = 1, i.e., classical MPC, cf. Figure 3. 

Examples 3.1 shows that choosing a larger value for to 
may significantly reduce the required prediction horizon 
length and, consequently, the numerical effort which grows 
rapidly with respect to the horizon N. 

Since the shape of the curve in Figure 2 is not a coinci- 
dence, i.e., aAr_|^jv/2j > C(N,m holds for to G {1, . . . , A^— 1}, 
cf. Griine et al. (2010), we may easily determine the small- 




Fig. 2. Visualization of the ai5^. -values with respect to 
the number of controls to be implemented m = 
1,..., [7V/2J} for prediction horizon N = 15. The 
dashed line indicates the reference value a = 0.275. 

est prediction horizon N G N>2 such that aN,[N/2\ > S 
ensures the desired performance specification. However, 




Fig. 3. Illustration of the a.,i-values with respect to the 
prediction horizon N G {15, . . . , 25} for m = 1, i.e., 
classical MFC. The dashed line indicates the reference 
value a = 0.275. 

note that these results exhibit a set-valued nature. Hence, 
the validity of the controllability condition (10) and the 
use of a prediction horizon length corresponding to the 
mentioned formula lead to an estimate which may be 
conservative at least in parts of the state space. This 
motivates the ckivclopment of algorithms which use the 
above calculated horizon length but close the control loop 
more often. Here, we make use of an extension of the m- 
stcp suboptimality estimate derived in Griinc and Pannek 
(2009) which is similar to Griine et al. (2010) but can be 
applied in a trajectory based setting. 

Proposition 3.2. Consider a G (0,1] to be fixed. If there 
exists a function Vjv : X — >• Mq satisfying 

mn — 1 

V^r(a;„) > Vjv(a;„+i) +a ^ £(a;„^(A;;a;„),ujv(fc;a;„)) 

(11) 

with m„ e {1, . . . , -/V - 1} for all n e No, then 

K^^(aro)< ^•Kx>(aro) (12) 
holds for all n e Nq. 
Proof: Reordering (11), we obtain 

m„ — 1 

a ^ i{Xu„{k\Xn),UN{k;Xn)) <VN{Xn) -VN(,Xn+l). 
fe=0 



Summing over n e No yields 

n rrii — l 

5^ 51 ^{XuN{hXi),UN{k\Xi)) 
i=0 fe=0 

< Vn{xq) - VAr(a;„+i) < VAr(a;o) 

and hence taking n to infinity implies the assertion. □ 

Rem,ark 3.3. Note that the assumptions of Proposition 3.2 
indeed imply the estimate 

Voo{Xn) < V^^ix,,) < ^VooiXn) (13) 

a 

for all n G No which can be proven analogously. 

Note that <S is built up during runtime of the algorithm 
and not known in advance. Hence, <S is always ordered. A 
corresponding implementation which aims at guaranteeing 
a fixed lower bound of the degree of suboptimality a takes 
the following form: 

Given state x := xo, list <S = (0), A'' e N>2, and a G (0, 1) 

(1) Set j = 0, compute un{-;x) and Vn{x). Do 

(a) Set j = J + 1, compute VAr(x^jj^, (j; a;)) 

(b) Compute maximal a to satisfy (11) 

(c) If a > a: Set to„ = j and goto 2 

(d) If J = N — 1: Print warning "Solution may 
diverge", set m„ = 1 and goto 2 

while a <a 

(2) For j = 1, . . . , m„ do 

Implement iJ.%{j — l;a;) := u^ij — ^',x) 

(3) Set S := {S, back(»S) + m„), x := x s (to„; x), goto 1 



Remark 3.4. Here, we adopted the programming notation 
back which allows for fast access to the last element of a 
list. 

Remark 3.5. If (11) is not satisfied for j < N — 1, a 
warning will be printed out because the performance 
bound a cannot be guaranteed. In order to cope with 
this issue, there exist remedies, e.g., one may increase the 
prediction horizon and repeat step 1. For sufficiently large 
N, this ensures the local validity of (11). Unfortunately, 
the proof of Proposition 3.2 cannot be applied in this 
context due to the prolongation of the horizon. Yet, it 
can be replaced by estimates from Pannek (2009) or 
Giselsson (2010) for varying prediction horizons to obtain a 
result similar to (12). Alternatively, one may continue with 
the algorithm. If there does not occur another warning, 
the algorithm guarantees the desired performance for a;„ 
instead of xq, i.e., from that point on, cf. Remark 3.3. 

In order to obtain a robust procedure, however, it is 
preferable that the control loop is closed as often as 
possible, i.e., m„ = 1 for all n e No. In the following 

section, we give resiilts and corresponding algorithms 
which given certain conditions allow us to reduce to„. 

4. ALGORITHMS AND STABILITY RESULTS 

In many applications a stable behavior of the closed loop is 
observed for m„ ~ 1 even if it cannot be guaranteed using 
the suboptimality estimate (6), cf. Example 2.2. Here, 
we present a methodology to close the control loop more 
often while maintaining stability. In particular, if to„ > 1 
is required in order to ensure a > a in (11) and, as a 



consequence, in step (1) of the proposed algorithm, we 
present conditions which aUow us to update the control 
for some j G {!,..., m„ — 1} to 

{^^f}..^ \ f'^N{k;Xn), k<j-l 
\uN{k ~ j;xu^^{j]x,i)), k>j 

insert the new updating time instant s(n) + j into the 
list iS, and still be able to guarantee the desired stability 
behavior of the closed loop. Summarizing, our aim consists 
of guaranteeing that the degree of suboptimality a is 
maintained, i.e., (11) is still satisfied, for the MPC law 
updated at a time instant j < to„ — 1 according to (14). 

Proposition 4-1- Let j:„ G X be given and inequality (11) 
hold for un{-; a;„), 5 > 0, and m„ e {2, . . . , A'' — 1}. If the 
inequality 

Vn{x {j;Xn)))- VN-j{Xuj^{j;Xn)) 



J-1 



< (1 - ") V ^i^UN {k; Xn), UN{k; Xn)) 



fc=0 



(15) 



-a ^ £( 

holds for some j <E {l,...,m„ — 1}, then the control 
sequence UN{-;Xn) can be replaced by (14) and subop- 
timality degree a is locally maintained. 

Proof: In order to show the assertion, wc need to show 
(11) for the modified control sequence (14). Reformulating 
(15) by shifting the running costs associated with the 
unchanged control to the left hand side of (15) we obtain 

(0;a;„)) 

— 1 

{k;Xn),UN{k;Xn)) 



k=0 



which is equivalent to 

VN{Xn) > VN{XiXf,{mn;Xn)) 

m„ — 1 



(16) 



-l-a 



fe=0 



^{Xun {k] Xn), Uf^ik] Xn)), 



i.e., the relaxed Lyapunov inequality (11) for the updated 
control zijv(s a^„). □ 

Example 4-2. Consider the illustrative Example 2.2 with 
horizon = 19 and a = 0.1. For this setting, we obtain 
a > a with m„ = 1 for all updating instances s„ G 5 
except for those three points indicated by (□) in Figure 1 
where we have rn„ = 2 once and m„ = 3 twice. Yet, for 
each of these points, inequality (15) holds with j = 1. 

• rn„ = 2: We apply Theorem 4.1 in order to update the 

control UN{'',Xn) '■= UAr(-;a;) according to (14). Thus, 
since the MPC law was updated at each sampling 
instant, we have obtained stability of the classical 
MPC closed loop. 

• m„ = 3: We utilize Theorem 4.1 in an iterative 
manner: since (15) holds for j — 1, wc proceed as 
in the case to„ = 2, update the control UN{-',Xn) ■— 
UN{-',Xn), and compute x^„(2;x„). Then, we com- 
pute Un{-; x^j^{2; a;„)) and check whether or not (15) 
holds for (14) with j = 2. Note that here, UN{-;Xn) 
in (14) is already an updated control sequence. Since 



(15) holds true, we obtain stability of the closed loop 
for the updated MPC control sequence. Moreover, we 
applied /Ltjv(x(s(n)), iJ,N{x{s{n) + l)) and ^j^{x{s{n) + 
2)) instead of iJ,ff{i;x{s{n))), i = Q, 1,2, i.e., we have 
again used classical MPC. 

Note that this iterative procedure can be extended in a 
similar manner for all m„ < N — 1. 

Integrating the results from Theorem 4.1 into the algo- 
rithm displayed after Proposition 3.2 can be done by 
changing Step (2). 



(2) For j = 1,...,TO„ do 

(a) Implement fiff{j — 1; x) := unU — h x) 

(b) Compute UN{-;Xuj^ij;x)) and l^Ar(a:ii„ (rn„; x)) 

(c) If j < m.n and condition (15) holds: 

Construct iiN{-',x) according to (14) 
Update ujv(-; x) := un{-; x) 



Note that VN-j{xu^{j; Xn)) in (15) is known in advance 
from VN{xn) due to the principle of optimality. Hence, 
only UN{-;Xu^{j;Xn)) and Vjv(a;u„(m„; x„)) need to be 
computed. Unfortunately, this result has to be checked for 
all j G {1, . . . , m„ — 2} using the comparison basis Viv(x„). 
As a consequence, we always have to keep the updating 
instant s{n) in mind, cf. Example 4.2. The following result 
allows us to weaken this restriction: 

Proposition 4- 3- Let x„ G X be given and inequality (11) 
hold for Un{-', Xn), 5 > 0, and m„ G {2, . . . , A'' — 1}. If the 
inequality 

{j;xn))) - l^jv(a;„„(m„;x„)) 

nin — l 

<a^ {£{xuj^{k;Xn),UN{k;Xn)) (17) 

k=j 

- £{xuN {k - j; Xuj,, (j; Xn)),UN{k - j; x„„ (j; x„)))) 

holds for some j G {1, . . . , m.„ — 1}, then the control 
sequence W]v(-;a;„) can be replaced by (14) and subop- 
timality degree a is locally maintained. 

Proof: Adding (17) to (11) leads to 
VAr(xfi„(m„;x„)) - Vn{ 

< -a^£{xuf,{k;Xn),UN{k;Xn)) 
fc=o 

{k;xn),UN{k;xn))- 

Updating the control sequence according to (14), we obtain 
(16) which completes the proof. □ 

Remark 4-4- We also like to point out that conditions 
(15) and (17) allow for a less fast decrease of energy 
along the closed loop, i.e., the case l^jv(2;u„(m„;x„)) > 
Viv(x„^(m„;x„)) is not excluded in general. 

As outlined in Example 4.2, Theorem 4.1 (and analo- 
gously Theorem 4.3) may be applied iteratively. Yet, in 
contrast to (15), VAr(x„) is not required in (17). Hence, 
if an appropriate update can be found for time instant 
j G {!,..., rrin — 1}, the loop can be closed and, as 
a consequence, Theorem 4.3 can be applied to the new 
initial value x<i(j; x„) with respect to the reduced number 



of control to be implemented ?n„ — j. To this end, the 
following modification of Step (2c) can be used: 



(2c) If j < nin and condition (17) holds: 

Construct un{-',x) according to (14) 
Update un{-;x) := un{-',x) 
Update the values VAr(a;„^(m„;a;)) and 
£{xuj^ {k; x), UN{k; a;)), A; = j, . . . , m„ - 1 



Combining the results from Proposition 3.2 and Theorems 
4.1, 4.3, we obtain our main stability result: 
Theorem 4-5. Let G X, cE > be given and apply 
one of the three proposed algorithms. Assume that (Ic) 
is satisfied for some j G {1, . . . , — 1} for each iterate. 
Then the closed loop trajectory satisfies the performance 
estimate (12) from Proposition 3.2. If, in addition, the 
conditions of Proposition 2.1(ii) hold, then Xjujv(-) behaves 
like a trajectory of an asymptotically stable system. 

Proof: The algorithms construct the set S. Since (Ic) is 

satisfied for some j G {1, . . . , N — 1} the assumptions of 



Proposition 3.2, i.e., inequality (11) for 



arc satisfied which implies (12). To conclude asymptoti- 
cally stable behavior of the closed loop, standard direct 
Lyapunov techniques can be applied. □ 
Exam,ple 4.6. Again consider Example 2.2 with A'^ = 19 
and a = 0.1. Since (17) is more conservative than (15), 
the closed loop control and accordingly the closed loop 



solution ; 



•) will not coincide in general. This also holds 



true in our example: While the case nin = 2 can be 
resolved similarly using the algorithm bascxl on (17), no 
update of the control ujv is performed for the two cases 
m„ = 3. Here, this does not cause any change in the list 
S. However, in general, «S is likely to be different for the 
presented algorithms. 

Since the proposed algorithms guarantee stability for the 
significantly smaller horizon A = 19, the computing time 
for Example 2.2 is reduced by almost 50% - despite the 
additional computational effort for Step (2). 

5. CONCLUSION 

Based on a detailed analysis of a relaxed Lyapunov in- 
equality, we have deduced conditions in order to close 
the gap between theoretical stability results and stable 
behavior of the MPC closed loop. Using this additional 
insight, we developed two MPC algorithms which check 
these conditions at runtime. Although the additional com- 
putational effort is not negligible, the proposed algorithms 
may allow us to shorten the prediction horizon used in the 
optimization. This reduces, in general, the complexity of 
the problem significantly which may imply a reduction of 
the overall computational costs. 

REFERENCES 

Alamir, M. and Bornard, G. (1995). Stability of a trun- 
cated infinite constrained receding horizon scheme: the 
general discrete nonlinear case. Automatica, 31(9), 
1353-1356. 

AUgower, F. and Zheng, A. (2000). Nonlinear model 
predictive control, volume 26 of Progress in Systems and 
Control Theory. Birkhauser Verlag, Basel. Papers from 
the workshop held in Ascona, June 2-6, 1998. 



Altmiiller, N., Griine, L., and Worthmann, K. (2010). 
Receding horizon optimal control for the wave equation. 
In Proceedings of the 49th IEEE Conference on Decision 
and Control, 3427-3432. Atlanta, Georgia. 

Badgwell, T. and Qin, S. (2003). A survey of industrial 
model predictive control technology. Control Engineer- 
ing Practice, 11, 733-764. 

Chen, H. and AUgower, F. (1999). Nonlinear model 
predictive control schemes with guaranteed stability. 
In Nonlinear Model Based Process Control, 465-494. 
Kluwer Academic Publishers, Dodrecht. 

Galaz, M., Ortega, R., Bazanella, A., and Stankovic, A. 
(2003). An energy-shaping approach to the design of 
excitation control of synchronous generators. Automat- 
ica J. IFAC, 39(1), 111-119. 

Giselsson, P. (2010). Adaptive Nonlinear Model Predictive 
Control with Suboptimality and Stability Guarantees. 
In Proceedings of the 49th Conference on Decision and 
Control, 3644-3649. Atlanta, GA. 

Graichen, K. and Kugi, A. (2010). Stability and incremen- 
tal improvement of suboptimal MPC without terminal 
constraints. IEEE Trans. Automat. Control. To appear. 

Grimm, G., Messina, M., Tuna, S., and Teel, A. (2005). 
Model predictive control: for want of a local control Lya- 
punov function, all is not lost. IEEE Trans. Automat. 
Control, 50(5), 546-558. 

Griine, L. (2009). Analysis and design of unconstrained 
nonlinear MPC schemes for finite and infinite dimen- 
sional systems. SIAM Journal on Control and Opti- 
mization, 48, 1206 1228. 

Griine, L. and Parmek, J. (2009). Practical NMPC sub- 
optimality estimates along trajectories. Sys. & Contr. 
Lett, 58(3), 161-168. 

Griine, L., Pannck, J., Seehafcr, M., and Worthmann, 
K. (2010). Analysis of unconstrained nonlinear MPC 
schemes with varying control horizon. SIAM J. Control 
Optim., 48(8), 4938-4962. 

Ito, K. and Kunisch, K. (2002). Receding horizon optimal 
control for infinite dimensional systems. ESAIM Control 
Optim. Calc. Var., 8, 741-760. 

Jadbabaie, A. and Hauser, J. (2005). On the stability of 
receding horizon control with a general terminal cost. 
IEEE Trans. Automat. Control, 50(5), 674-678. 

Keerthi, S. and Gilbert, E. (1988). Optimal infinite- 
horizon feedback laws for a general class of constrained 
discrete-time systems: stability and moving-horizon ap- 
proximations. J. Optim. Theory AppL, 57(2), 265-293. 

Magni, L., Raimondo, D., and AUgower, F. (2009). Non- 
linear model predictive control. Towards new challenging 
applications. Springer Berlin. 

Magni, L. and Scattolini, R. (2007). Robustness and 
robust design of MPC for nonlinear discrete-time sys- 
tems. In R. Findeisen, F. AUgower, and L.T. Biegler 
(eds.). Assessment and future directions of nonlinear 
model predictive control, volume 358 of Lecture Notes 
in Control and Inform. Sci., 239-254. Springer, Berlin. 

Pannck, J. (2009). Receding Horizon Control: A 
Suboptimality based Approach. PhD Thesis in Math- 
ematics, University of Bayreuth, Germany. 

Rawlings, J.B. and Mayne, D.Q. (2009). Model Predictive 
Control: Theory and Design. Nob Hill Publishing, 
Madison. 



